Parity- Affected Superconductivity in Ultrasmall Metallic Grains 
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We investigate the breakdown of BCS superconductivity in ultrasmall metallic grains as a function 
of particle size (characterized by the mean spacing d between discrete electronic eigenstates), and 
the parity (P = even/odd) of the number of electrons on the island. Assuming equally spaced 
levels, we solve the parity-dependent BCS gap equation for the order parameter Ap(d, T). Both the 
T = critical level spacing d Cj p and the critical temperature T Ct p(d) at which Ap = are parity 
dependent, and both are so much smaller in the odd than the even case that these differences should 
be measurable in current experiments. 



PACS numbers: 74.20.Fg, 74.80.Fp, 74.80.Bj 

The study of the properties of ultrasmall metallic par- 
ticles has witnessed a dramatic development during the 
last year: using an ingenious new fabrication technique, 
Black, Ralph and Tinkham (BRT) 111 have constructed 
a single-electron transistor (SET) whose island, a single 
nm-scale Al grain, is more than four orders of magnitude 
smaller in volume (estimated radii between r ~ 2.5 and 
13 nm) than that of conventional SETs. Thus a new 
energy scale, the average level spacing d = 1/N(ep) be- 
tween discrete electronic levels, enters the problem: Both 
the free-electron estimate of d ~ 2-k 2 Ti 2 / {mkp V) and di- 
rect observation (discrete steps in the I-V curve) give 
values of d ranging from 0.02 to 0.3 meV, the latter be- 
ing much larger than the smallest accessible temperatures 
(~ 30 mK) and on the order of the bulk superconducting 
gap (A b = 0.18 meV for Al). 

The eigenenergies of the larger grains (r > 5 nm) stud- 
ied by BRT revealed the existence of an excitation gap 
O > d which is driven continuously to zero by an ap- 
plied magnetic field, and striking gap-dependent parity 
effects, i.e. differences between islands with an even or 
odd [P = e/o] number of electrons. BRT very convinc- 
ingly interpreted these phenomena as evidence for su- 
perconductivity. On the other hand, in smaller particles 
(r < 5 nm) it was not possible to detect an Jf-dependent 
gap or other evidence for superconductivity. 

These experiments invite reconsideration of an old but 
fundamental question: what is the lower size limit for 
the existence of superconductivity in small grains? An- 
derson addressed this question already in 1959 Q and 
argued that "superconductivity would no longer be pos- 
sible" if the level spacing d becomes larger than the bulk 
gap Ab, for reasons explained below. This answer - al- 
though in general correct - is not yet quite complete, 
since it does not address parity effects. Even in "large" 
superconducting islands (with d <C At) experiments || 
have demonstrated the dramatic impact of parity on I-V 
characteristics; moreover theory predicts an even- 
odd difference for the superconducting order parameter 
itself of A e -A = d/2 at T = 0. Though the latter dif- 



ference is immeasurably small in "large" islands, it should 
certainly become significant in ultrasmall grains. More- 
over, since the crossover temperature at which parity ef- 
fects become observable || , namely T cr — Ab/ In N e f f 
(where in the d < A b limit N eff = ^/8irTA b /d), be- 
comes of order Ab when d ~ Ab, parity effects should sur- 
vive to temperatures as high as T c itself. Hence T c p(d) 
as function of d should be parity-dependent too. 

In this Letter we address these issues by studying par- 
ity effects in the order parameter Ap(d, T) for general d. 
In particular, we calculate Ap(d, 0) and T Cy p{d) by solv- 
ing the BCS gap equation (derived using parity-projected 
mean- field theory j|]|]) at T — and Ap = 0, respectively, 
for the case of equally-spaced single-particle levels. We 
find T C)0 (d)/T C;e (d) < 1 and a remarkably small ratio 
of critical level spacings d Ct0 /d c ,e — 1/4 at T = 0. Our 
results are completely compatible with BRT's observa- 
tions. Moreover, the predicted parity differences should 
be directly observable in their latest experiments which 
have variable gate voltage, allowing them to change the 
number parity of a given grain at will. 

The model: In BRT's experiments, the charging energy 
Eq = e 2 /2C 'total of an ultrasmall grain is by far the largest 
energy scale in the problem (with Ec — 4 meV ^> Ab), 
so that fluctuations in particle number are strongly sup- 
pressed. Therefore, in this Letter we consider a com- 
pletely isolated grain, which should be described using a 
canonical ensemble with a prescribed number of electrons 
n = 2m + p, where p = (0, 1) for P = (e, o) (the labels p, 
P and also n will be used interchangeably as parity la- 
bels below). We adopt a model Hamiltonian having the 
standard reduced BCS form: 
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(1) 



Here cj CT creates an electron in the particle-in-a-box-like, 
independent-electron state \ja), where the states \j+) 
and \j—) are degenerate, time- reversed partners whose 
energies are considered as a given set of phenomeno- 
logical parameters. The integer j is a discrete quantum 
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number. For a given n = 2m + p, we take j = to de- 
scribe the first energy level whose occupation in the T = 
Fermi sea is not 2 but p, so that j = —m, . . . oo. Finally, 
the dimensionless coupling constant A -1 = In =^ is re- 
garded as a phenomenological parameter determined by 
the value A = A(0, 0) of the effective gap (measured at 
d<C A) and some cut-off frequency lu c . 

Pair-mixing: At this point it seems appropriate to 
briefly address the question of what is meant by the "ex- 
istence of superconductivity" in ultrasmall grains. It de- 
serves special attention, firstly because the BCS order pa- 
rameter A, if defined as usual by Ad£^ - (c;_Cj+), seems 
to be zero in a canonical ensemble, and secondly because 
most of the standard criteria, e.g. a gap followed by a 
continuous excitation spectrum, zero resistivity and the 
Meissner effect, are not applicable here. 

Now, the microscopic reason for all of these (large- 
sample) phenomena is of course the existence of a pair- 
correlated ground state. The essence of its correlations 
is what we shall call pair-mixing across ef, namely the 
partial population (with amplitude Vj > 0) of some time- 
reversed pairs of states above ef {j > 0) 
by partially depopulating (with amplitude uj > 0) some 
pairs of states below £p (j <0). This creates phase space 
for pair scattering (which is Pauli-blocked in the normal 
ground state) and hence allows the BCS interaction to 
lower the ground state energy. 

Although BCS showed that a brilliantly simple way of 
calculating the Uj and Vj is to use grand-canonical meth- 
ods, pair-mixing of course can and does also occur in a 
fixed-n system. Indeed, this pair-mixing can readily be 
characterized by a "generalized" order parameter that is 
equal to the conventional ^dJ2j( c j- c j+) m BCS's grand- 
canonical mean- field treatment, but (in contrast to the 
latter expression) is meaningful in a fixed-n system too, 
namely Arf^j UjVj. An experimental signature of this 
pair-mixing is the energy cost needed to add or remove 
single electrons that perturb these correlations (i.e. that 
"break pairs"). Since BRT quite unambiguously mea- 
sured such energy costs in their larger grains, it seems 
reasonable to regard these as "superconducting" , in the 
sense of having a pair- correlated ground state that mea- 
surably exhibits pair-mixing. 

The notion of pair-mixing also provides a simple way 
to understand why superconductivity ceases to exist in 
sufficiently small samples. If the level spacing becomes 
sufficiently large (d ~ A) , pair-mixing costs a prohibit- 
ive amount of kinetic energy and hence ceases to occur. 
The task at hand is to describe this breakdown (semi)- 
quantitatively, while keeping track of parity effects. 

Canonical and Parity Projection: Since in practice it 
is so much harder to calculate Uj, Vj canonically than 
grand-canonically, and since the latter approach has been 
remarkably successful even in small systems such as nu- 
clei with as little asn~ 100 nucleons B , purely canonical 
descriptions (pioneered by Schafroth M) are seldom used. 



An alternative approach |^,|| is to employ an auxiliary 
parity-projected grand-canonical partition function, 

Z${ l ,)=Tr G \[l±{-f]e-^-^ = e-^ZM , (2) 

(Tr G denotes a grand-canonical trace), from which the 
desired fixed-n partition function Z n can in principle be 
exactly projected by integration: 

Z n = f ^e- im Z${iu/p) . (3) 

Since in practice, though, it is hard to perform the inte- 
gration exactly, we approximate the integral by its saddle 
point value, Z n ~e~P fJ " nTl Zp(fi n ), where fi n is fixed by 

n = -^0«(/i)| M=Atn [={N)p] (4) 

(Here ( )p is taken in the parity-projected grand- 
canonical ensemble of Zp.) This equation, the brack- 
eted part of which is the parity-projected version of a 
standard grand-canonical identity, illustrates the elemen- 
tary fact that the saddle point approximation produces 
nothing but the grand-canonical description we had set 
out to improve upon. Nevertheless, the merits of the 
above approach are that (i) corrections can be calculated 
systematically by including Gaussian contributions [gov- 
erned by ((TV 2 ) - (AO 2 )" 1 / 2 ] to the integral (§) (at T = 
they are ~ (d/A) 1 / 2 W); (ii) it illustrates that the par- 
ity projection of Eq. (0), which is essential for extracting 
e/o differences, can be done exactly even when the fixed- 
n projection cannot; (iii) it clarifies that in a canonical 
ensemble \x n is simply the saddle-point value of an inte- 
gration parameter, which, however, has to be determined 
with special care in ultrasmall grains, for which d is large. 

Mean-Field Approximation: We evaluate Zp using a 
"naive mean-field approach" (our method is equivalent 
to that used in ||): Make the replacement 

C j- C 3 + i C J- C 3+ - ( C i- C j+)P} + ( C 3- C 3 + )P ( 5 ) 

in H — n n N, neglect terms quadratic in the fluctua- 
tions represented by { } and diagonalize, using J n j<7 = 
UnjCnja — Gv njc\^_ a . One obtains the usual results 
H-^nN ~ C n + J2 jlT Kj^lj^nja, where E njr7 = [e^ + 
Ap] 1/2 , e nj =e^ - ^ and vf l} =|(1 - e nj /E nj ). More- 
over, since the parity of electron number and quasiparti- 
cle number are always the same, Eq. (|^) can be rewrit- 
ten H using quasiparticle-parity projection, Zp(fi n ) = 
\{ZX±Z a ), where 

Z°(fi n )=e-^l[(l±e-^). (6) 

3<r 

The usual mean-field self-consistency condition Ap = 
AdJ\ - (cj_Cj+)p takes the form 
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njcr 7 



(7) 



where f nja = (l nj(7 lnja)p is given by ffl 
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Here we defined / 



The above description thus involves the usual BCS quasi- 
particles, but their number parity is restricted to be p. 

Corrections to the naive mean-field expressions for 
E n j a and Ap [as measured, e.g., by fluctuations in Ap, 
or by (AdS^f^p)" 1 / 2 ] can be shown to be of order Xd /A. 
Janko et al. Q] calculated the leading one of these correc- 
tions (their 1 /N e f f terms) by using a more careful mean- 
field approach which incorporates contributions from the 
{ } 2 terms neglected above. However, for Xd/ A suffi- 
ciently large that these terms matter, other non-mean- 
field corrections (~ (Xd/ A) 1 / 2 ) matter too. Neverthe- 
less, since our goal is merely to find a criterion for when 
pair-mixing correlations cease to exist, and one such cri- 
terion is to establish when the mean-field gap equation 
no longer has a solution, we shall neglect all corrections 
going beyond naive mean-field theory 

Determination of fi n : Following Q , let us henceforth 
consider the case of equal level spacing, e® = j d 
(which seems reasonable for large n, due to level repul- 
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sion). Using (N) P = J2ja { v l 3 + ( u: 
symmetrizing the j-summation interval by dropping the 
negligibly small j>m— <5p e terms, Eq. (Eh gives 
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Here we have added and subtracted sign(e„j) (defin- 
ing sign(O) = and 6(0) = |). By inspection, this implies 
that fj, n = £q — ^dSp te , because then e n j = d(j + \5p : e) , 
which ensures firstly that 29(e n o) = p, and secondly that 
the summands of the second and third terms are anti- 
symmetric, so that these terms vanish. (More generally, 
they vanish also for non-equally-spaced levels as long as 
the levels are distributed roughly symmetrically about 
Eq, in which case /x„ ~ £o<5p,o + ^(^o — £ °-i)°~P,e) Thus 
we have proven the seemingly obvious: in the language 
of the T = normal Fermi sea, /i„ lies exactly half-way 
between the last filled and first empty levels if P = even, 
and exactly on the singly occupied level if P— odd. 

We are now ready to study the gap equation (Q). 

Gap Equation at T = 0: The quasiparticle occupation 
function reduces to f n j a — ^<5jo^p.o at T = 0, as intu- 
itively expected, because then the even or odd systems 
have exactly zero or one quasiparticle, the latter in the 
lowest quasiparticle state, namely j — 0. This e/o differ- 
ence has a strong impact on the T = gap equation: in 



the odd case, the j = level, for which E~- is largest, 
is absent, reflecting the fact that the odd quasiparticle 
in the j = state obstructs pair scattering involving this 
state. To compensate this missing term, A G must there- 
fore become significantly smaller than A e as soon as d 
is large enough that a single term becomes significant 
relative to the complete sum. 

To quantify this statement, it is convenient to rewrite 
Eq. (0) at T = as follows: Writing E~j = J dtu/ir (E nj + 
lu 2 )^ 1 , transferring the cut-off w c from to Jdu>, and 
performing the j-sum (by contour integration) gives 
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where Ep u = (uj 2 + Ap) 1 / 2 . Since, amusingly, for P =e 
Eq. (||) is identical in form (with d — > 2nT) to the well- 
known gap equation for the T-dependence of the bulk 
gap (curve A in Fig. 0), we have A e (d, 0) = Ap(0, d/2ir). 
Thus, for d/A <C 1 the even gap has the standard form 



A e (d,0) = A(l- 



-2nA/d\ - n cori |; ra gt, one eas- 



ily finds from Eq. ([|) that the odd gap drops linearly, 
A Q (d, 0) = A - d/2, in agreement with ||J|. 

The full solutions of Eq. (0) for Ap(d, 0), obtained nu- 
merically, are shown as curves B and C in Fig. @. The 
critical values d c .p at which Ap(<i c .p, 0) = can be found 
analytically by setting Ap = T = in Eq. (0): 

' * ' " ' \e~<~ 0.890. (10) 
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Critical temperature: Although ultrasmall grains can- 
not undergo a sharp thermodynamic phase transition 
(this would require n — > oo), the quantity T Ct p(d), de- 
fined simply as the solution to the Ap — > limit of 
El- (0); is another measure of how rapidly pair-mixing 
correlations break down as function of level spacing. Our 
numerical results for T Cj p(d) JIo[| , shown as curves D and 
E of Fig. |l| for P — e/o, have the expected limits at d = 
and G? c ,p, but behave unexpectedly in between: 

Even: In the even case, T Cje (d) is non-monotonic, in- 
creasing slowly from its initial value T c e (0) = Ae 7 /7r = 
0.567A to a maximum given by T c , e (2.60A) = 0.736A > 
T Cie (0). It then drops to zero very rapidly as T c ^ e (d) ~ 
dj log[4d Ci e/ (d c ,e — d)], which follows from setting A e = 

in A e (d,T) ~ d[2(l-d/d c , e -4e- d / T )/(7C(3))] 1/2 , an 
analytical result that holds for the regime T, A e (T) <C d. 
The intuitive reason for the initial increase in T c ^ e (d) 
is that the difference between the actual and the usual 
quasiparticle occupation functions is fnja~fnja < for an 
even grain (becoming significant when d ~ A), reflecting 
the fact that exciting quasiparticles two at a time is more 
difficult than one at a time. Therefore the quasiparticle- 
induced break-down of superconductivity with increasing 
T will set in at slightly higher T if d ~ A. 

Odd: In the odd case, T CjD (d) "overshoots" a bit near 
d = d c<0 , or alternatively, the critical level spacing d Cj0 (T) 
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is non-monotonic as a function of increasing T, reaching 
a maximum value given by d CiO (0.387A) = 1.05A before 
beginning to decrease. The intuitive reason for this is 
that for < A < T, d the odd j = function /„ 0ct (T) 
becomes somewhat smaller than its T = value of i, 
because with increasing T some of the probability for 
finding a quasiparticle in state j "leaks" from j = to 
higher states with j ^ 0, for which E~j < E~} in Eq. (0). 
Thus, the dramatic blocking-of-pair-scattering effect of 
the odd quasiparticle becomes slightly less dramatic as 
T is increased, so that d c ^ a increases slightly. 

Discussion: An important general feature of our re- 
sults is that level discreteness always reduces Ap{d, 0) to 
be < A (thus contradicting Ref. || , which was convinc- 
ingly criticised in Ref. §). However, BRT's experiment 
found an effective gap A that is larger by a factor of 1.5 to 
2 than its bulk value A&. Following the argumentation 
of H for thin films, wc can attribute this to presumed 
changes in the phonon spectrum in small samples, which 
can be modeled by using a constant value of A larger (by 
a few percent) than the usual bulk value \t- 

The rather rapid drop of Ap (d) , once it happens, could 
be the reason why BRT see a well-developed gap A 
even for d ~ A, but don't see any for their smallest 
grains. The most important conclusion of this paper, 
though, is summarized by Fig. |l| and Eq. (]To)) : there is 
a large regime in which A D = while A e is still ~ A, 
in other words, pair-mixing correlations vanish signifi- 
cantly sooner for odd than even grains as their size is 
reduced. Moreover, the surprising largeness of the ratio 
d c ,e/d c .o — 4 opens the exciting possibility of studying 
grains with d c ^ < d < d c ,ei which should have A D = 
while A e is still ~ A. Since this leaves tell-tale traces 
in the I-V-V g characteristics jL0|, BRT should be able 
to test this prediction directly, since they can change the 
electron parity of a given grain in a controlled way by 
tuning the gate voltage of their SET, and hence study 
the same grain in both its even and odd states. 

Of course, a nagging question remains: How robust are 
our mean-field-based results? Will not fluctuations in A, 
governed by d/A and hence large in ultrasmall grains, 
wash out the predicted parity differences in A e / Q ? For in- 
stance, it is doubtful that the unexpected non-monotonic 
subtleties of T Ci p(<f), though intuitively plausible, have 
physical significance, since they fall in the Ap ~ regime 
where fluctuations are extremely large. 

However, while a detailed analysis of the fluctuations is 
beyond the scope of this Letter, we believe that, at least 
in the (experimentally accessible) regime of T/d ~ 0, our 
main result is indeed robust: Firstly, it is known that 
at T = 0, the leading d/A fluctuation corrections to the 
free energy difference fi^f — are parity-independent 
||,[|. Secondly, since as d approaches the point where 
A D vanishes, we have d/A e < 1 <C d/A a , fluctuations 
become important much later in the even than the odd 



case. Thirdly, even for systems much smaller than ul- 
trasmall grains (that have n ~ 10 4 ), namely shell model 
nuclei (with n ~ 100), the T = BCS-description of 
pairing interactions has been remarkably successful 0, 
despite the presence of large fluctuations. 

The perhaps most compelling support for the reality 
of the predicted effect comes from the following back-of- 
the-envelope argument: the ground state energy differ- 
ence between a BCS superconductor with parity P and 
the normal Fermi sea is — A P /(2d) + A Sp. . While for 
P — e it is always negative, for P = o it vanishes when 
d> 4 A, illustrating that the odd electron indeed dramat- 
ically disrupts pair-mixing correlations when d~A. 

In conclusion, we have investigated the influence of 
parity on the superconducting mean-field order parame- 
ter in ultrasmall grains. We have found that as a function 
of decreasing grain size, superconductivity breaks down 
in an odd grain significantly earlier than in an even grain, 
which should be observable in present experiments. 

It is a pleasure to thank BRT for showing us their 
preliminary results, and to acknowledge discussions with 
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